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Abstract 



ff) ■ We consider the Adler D-function, which is defined in the spacelike region, 

and the e + e~ annihilation ratio smeared, according to the Poggio, Quinn, and 
\ Weinberg method, which is determined for timelike argument. We show that 

the method of the nonperturbative a-expansion allows one to describe these 
Euclidean and Minkowskian characteristics of the process of e + e~ annihilation 
down to the lowest energy scales. 

1. In quantum chromodynamics, it is important to determine the "simplest" objects which 
allow one to check direct consequences of the theory without using model assumptions in an 
essential manner. Comparison of theoretical results for these objects with experimental data 
allows us to justify transparently the validity of basic statements of the theory, and make some 
conclusions about completeness and efficiency of the theoretical methods used. Some single- 
argument functions which have a straightforward connection with experimentally measured 
quantities can play the role of these objects. A theoretical description of inclusive processes 
can be expressed in terms of functions of this sort. Let us mention among them moments 
M n (Q 2 ) of the structure functions in inelastic lepton-hadron scattering and the hadronic 
correlator Il(s) (or the corresponding Adler D-function), which appear in the processes of 
e + e~ annihilation into hadrons or the inclusive decay of the r lepton. 

In this paper we consider such objects for the process of e + e~ annihilation into hadrons. 
Here, the point is that the experimentally measured ratio of hadronic and leptonic cross- 
sections, R e + e - (s), is not suitable at the present stage of development of the theory to provide 
a description independent of model considerations. However, it is possible to construct simpler 
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objects than the _R e + e --ratio. These are the Adler function D(Q 2 ) IJ, which is defined in the 
Euclidean region, and the quantity Ra(s) constructed by the Poggio, Quinn, and Weinberg 
"smearing" method || and defined in the Minkowskian region. 

The theoretical method, which we will use here, is the nonperturbative expansion tech- 
nique suggested in Refs. 0. This approach is based on the idea of variational perturbation 
theory (VPT) H, which in the case of QCD leads to a new small expansion parameter. To 
compare the results with experiment we will use the new "experimental data" for the Adler 
function that has been recently obtained in Ref. |J and the smeared "experimental curve" 
corresponding to the i?-ratio taken from Refs. 

2. We use the method of constructing the so-called floating or variational series in quantum 
theory. Within this approach, a certain variational procedure is combined with the possibility 
of calculating corrections to the principal contribution which allows the possibility of probing 
the validity of the leading contribution and the region of applicability of the results obtained. 
At present, this idea finds many applications in the development of various approaches, which 
should enable us to go beyond perturbation theory. Among these are the Gaussian effective 
potential method ||, the optimized 5-expansion ||, and the VPT approach ||. 

We will apply a nonperturbative QCD expansion based on a new small expansion param- 
eter ||. Within this method, a quantity under consideration can be represented in the form 
of a series, which is different from the conventional perturbative expansion and can be used to 
go beyond the weak-coupling regime. This allows one to deal with considerably lower energies 
than in the case of perturbation theory. 

The new expansion parameter a is connected with the initial coupling constant g by the 
relation 

g 2 ots_ _ 1 a 2 m 
" (4vr) 2 " 4tt " C (1-a) 3 ' [ ' 

where C is a positive constant. As follows from (|l|), for any value of the coupling constant g, 
the expansion parameter a obeys the inequality 

< a < 1 . (2) 

While remaining within the range of applicability of the a-expansion, one can deal with low- 
energy processes where g is no longer small. 

The positive parameter C plays the role of an auxiliary parameter of a variational type, 
which is associated with the use of a floating series. The original quantity, which is approx- 
imated by this expansion, does not depend on the parameter C; however, any finite approx- 
imation depends on it due to the truncation of the series. Here we will fix this parameter 
using some further information, coming from the potential approach to meson spectroscopy. 
In the framework of this approach consider the following approximations to the renormal- 
ization group /^-function, the functions (3^ and which are obtained if one takes into 
consideration the terms 0(a 3 ) and 0(a 5 ) in the corresponding renormalization constant Z\. 
As has been shown in Ref. C is determined by requiring that -fr k \\)/\ tends to 1 for 
sufficiently large A, which gives C3 = 4.1 and C5 = 21.5. The increase of with the order of 
the expansion is explained by the necessity to compensate for the higher order contributions. 
A similar phenomenon takes place also in zero- and one-dimensional models. The behavior 
of the functions — /3^(A)/A gives evidence for the convergence of the results, in accordance 
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with the phenomenon of induced convergence. [] The behavior of the /3-function at large value 
of the coupling constant, —f3^ k \\)/\ ~ 1, corresponds to the infrared singularity of the run- 
ning coupling: a s (Q 2 ) ~ Q~ 2 at small Q 2 . In the potential quark model this Q 2 behavior is 
associated with the linear growth of the quark-antiquark potential. 

The renormalization group /^-function of the expansion parameter a is 



2 da 2(3 1 



where (3q — 11 — 2//3 is the one-loop coefficient of the /3-function in the usual perturbative 
expansion, and / is the number of active quarks, has a zero at a = 1 that demonstrates the 
existence of the infrared fixed point of the expansion parameter and its freezing-like behavior in 
the infrared region. By finding the renormalization constants in the massless renormalization 
scheme with an accuracy 0(a 3 ), we find for the function F(a) 

„,»s, 2 6 , 18 1 624 . . 5184 / 9 \ 

F< 8 >(o = 48 1na - - + —In 1 - a + — - In l + -o . 4 

cr a 11 1 — a 121 121 V 2 / 

By solving the renormalization group equation (||) we find the momentum dependence of 
the running expansion parameter a(Q 2 ) as a solution of the following transcendental equation 



laS = ^[F(o)-F(oo)]. (5) 



For any values of Q 2 , this equation has a unique solution a = a(Q 2 ) in the interval between 
and 1. 

By working at 0(a 5 ) we obtain a more complicated result 

F^(a) = 5{5 l w) I>J(aA) (6) 

with B = Pi/(2CPo), where the two-loop coefficient /3i = 102 — 3//3, and 

2 4 _ 12 9 4 + 126 + 216 2 

(G ' } ~ ~ ~ d> ~ (l-a)(l-6) + P a 

3 °- 21& , „ \ (2 + &) 2 w 

+ Ff ln(1 " fl) "^rf ln(fl ^ ) (7) 

with 



(bi - bj)(bi - b k ) 



Here indices {ijk} are {123} and cyclic permutations. The values of foj are the roots of the 
equation ip(bi) = 0, where the function i/j(o) is related to the /3-function and is 



It has been observed empirically [ 10 1 that the results seem to converge if the variational param- 
eter is chosen, in each order, according to some variational principle. This induced-convergence 
phenomenon is also discussed in Ref. 
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ij,{a) = 1 + + 2(6 + a)a 2 + 5(5 + 3 B)a 3 . (9) 

3. The cross-section for the process of e + e~ annihilation into hadrons, or its ratio to the 
leptonic cross-section, R(s), is a physically measured quantity, defined for timelike momentum 
transfer - the physical region of the process. These quantities have a resonance structure 
that is difficult to describe without model considerations. Moreover, the basic method of 
performing calculations in quantum field theory, perturbation theory, becomes ill-defined due 
to the so-called threshold singularities. Both these problems can, in principle, be avoided if 
one considers a "smeared" quantity. In the paper , Poggio, Quinn, and Weinberg suggested 
that instead of using the ratio 

i?( s ) = I[n( S + ie)-n( S -ie)] , (10) 

defined through the hadronic correlation function n(g 2 ) taken near the cut, which has a 
problem with threshold singularities, the following smeared quantity be used: 

R A (s) = i [II(s + iA) - U(s - iA)] (11) 

with some finite A. It has been argued in Ref. that for appropriate values of this smearing 
parameter (A is of order of a few GeV 2 ) it is possible to compare theoretical predictions with 
smeared experimental data. To get these data one can use the dispersion relation for the 
hadronic correlator Tl(s) together with Eq. (|ll|) to give 

Ra(s) = - T ds' R ^ MA2 . (12) 
7i Jo [s — s'Y + A z 

The corresponding "experimental" curves of R&(s) have been given in Refs. |||7]]. We will 
use these data to compare our results with experiment. 

However, a straightforward usage of conventional perturbation theory is not still possible. 
Indeed, by parametrizing, as usual, the QCD contribution to the function R(s) in Eq. flT2] ) 
by the perturbative running coupling, which has unphysical singularities, one encounters 
difficulty with the definition of the integral on the right-hand side. Moreover, the usual method 
of the renormalization group gives a Q 2 -evolution law of the running coupling in the Euclidean 
region and there is the question of how to parametrize, in terms of the same scale parameter 
A, a quantity, for example R(s), defined for timelike momentum transfer. Here an important 
role is played by the analytic properties of the running coupling. Within the nonperturbative 
a-expansion it is possible to maintain such properties and self-consistently determine the 



effective coupling in the Minkowskian region fL2| . Note that the so-called analytic approach 



to QCD [TJ| also leads to a well-defined procedure of analytic continuation from the spacelike 
to the timelike domain |H,15fl. By using the analytic approach the characteristics of the e + e~ 



annihilation process have been analyzed in Refs. ]T6] , ^7| . 

Another function which characterizes the process of e + e~ annihilation into hadrons is the 
Adler function 
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Being denned outside of the resonance region, the D-function is the more useful object to test 
QCD both in the ultraviolet perturbative region of large Q 2 and in the nonperturbative region 
of small Q 2 . Recently a new "experimental" curve for this function has been obtained ||. 
This curve is a smooth function of Q 2 without any traces of resonance structure, which makes 
it useful for comparison with reasonable theoretical descriptions. We will now consider D(Q 2 ) 
and Ra(s) in the framework of the nonperturbative a-expansion. 
In the massless case, let us represent these functions in the form 

D(Q 2 )=3Y,q 2 f [l + d \ cS (Q 2 )] (14) 
/ 

and 

R(s) = 3j2q 2 f [l+ro\f(s)} . (15) 
/ 

Here qj are the quark charges, do and r$ are the first perturbative coefficients which are 
renormalization scheme independent, being do = ro = 4. For the D-function defined in the 
spacelike region we can write down the effective coupling in the form of the a-expansion which 
at the 0(a 5 ) level has the form 

,nff 1 9 3 , f 6 1 di\ , ( 10 6 di\ ^ ,„„ N 

A .* = _ a * + _ a 3 + (_ + __j a < + (_ + __j a5 , (16) 

where d\ is the next coefficient of the perturbative representation for the D-function which 
we will take in the MS renormalization scheme ]19| 



d x = ~ [365 - 22/ - 8C(3) (33 - 2/) ] ~ 16 (1.986 - 0.115/) . (17) 

Here ((n) is the Riemann ^-function, £(3) ~ 1.202. 

Eq. (|15|) serves to define the effective coupling in the timelike domain or as we will say in 
the s-channel, which is reflected in the subscript s. For self-consistency of Eqs. (|i~4D and fli~5|) 
with the dispersion representation in Eq. (p~3| ) it is important to maintain the correct analytic 



properties of the initial effective coupling A (Q 2 ) |l^JT^JT5|| . In this case one finds 



A eff (Q 2 ) =Q 2 f 7-^71 Af(.) (« 
Jo (s + Q 2 ) 



and the corresponding inverse relation 



z =s 



The s-channel running coupling can be written in the form 



A5°(a) = (20) 

Z7T1 I po 



where a± obey the equation 
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F{a ± ) = F(a ) + 



c 



In 



Ql 



± i 7T 



At the level 0(a 3 ) the function <j>(a) has the form 
0( 3 ) (a) = -4 In a 



72 1 318, M . 256, / 9 

1 ln(l - a) H In 1 + -a 

111 -a 121 v ; 363 V 2 



(21) 



(22) 



Similarly, a more complicated expression for the 0(a 5 ) level which we will use can be derived. 
Note here that, as it has been recently argued from general priciples, the behavior of the 



effective couplings in the spacelike and the timelike regions cannot be symmetrical |T8 



To incorporate quark mass effects, we use for the cross-section ratio an approximate ex- 
pression 

R(s) = 3]TgJ©(s- Am}) T(v f ) [1 + g(v f )r f (s)} , (23) 
/ 

where the sum is performed over quark flavor and 

v f : 



'1- 



4m 



/ 



T(v) 



v(3-v 2 ) 



9{v) = — 



7T 



3 + v 



7T 



3 

47T 



_2v 4 V2 

The quantity r/(s) is defined by the s-channel effective coupling A eff ( 



(24) 

The corresponding 



-D-function can be found from Eq. (|13|). 

For massless, MS-like, renormalization schemes one has to consider some matching proce- 
dure. To this end, one usually adopts a procedure for matching the running coupling in the 
Euclidean region f20j by requiring that the running couplings, corresponding to (/ — 1) and 
/ number of fermions, should coincide with each other at some matching point Q = £ra/. 
For the matching parameter £ one usually takes 1 < £ < 2. Obviously, in this case the 
derivative of the running coupling will not be a continuous function and the correct analytic 
properties of the D-function will be violated. Instead we will apply the matching procedure 
in the s-channel, where, at least in the leading order, the number of active quarks is directly 
connected with the energy y/s and discontinuously changes at the threshold s = 4m^.| The 
effective charge in the Euclidean region, restored by the dispersion relation (|18|) , will have the 
correct analytic properties and will "know," in a way similar to massive MOM renormalization 
schemes, about all physical thresholds. 

To perform this matching procedure, one can require that the s-channel running coupling 
and its derivative be continuous functions in the vicinity of the threshold. This requirement 
leads to a system of equations, which we write down for a simple 0(a 3 ) case: 



/-i 



1% 



Im 



Po 



(25) 



Cf-i 



-Im 



3a+ 1 



— Im 



1 + 3a 



2 In the framework of the analytic approach, where there is also a well-defined procedure of analytic 
continuation, it is possible to use the same s-channel matching method fl5||. 
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The function 0(a) is denned by Eq. fl22|) and a + obeys to Eq. (|21[) . Therefore, Eqs. (p5[) allows 
one to establish relations between the parameters C/_i , ag" 1 and C/ , a{. 

Our results are presented in Figs. [l|-44|. In Figs. |l] and ||] we plot the smeared functions 
Ra(s) for A = 1 GeV 2 and A = 3 GeV 2 respectively. The solid line is the VPT next-to-lea- 
ding order (NLO) result which was normalized at the r lepton scale. To this end, we use the 
method of description of the r decay suggested in Refs. pl|j2^] and take here the following 
experimental average value R T = 3.642 ± 0.021 |23| as input. To demonstrate the fact of 
stability we show in Fig. ||] the leading-order (LO) VPT result. In these figures we also plot, 
as dashed curves, results obtained by using the principle of minimal sensitivity (PMS) to 
optimize the third order of the perturbative expansion^] and the smeared experimental data, 
as dotted lines, taken from Ref. |J. For A = 1 GeV 2 the resonances in the e + e~ annihilation 
cross-section are not smeared enough, since there are some peaks in the region of J / ip meson 
family. With increase in the value of A the resonance structure is smoothed and finally 
disappears. For a wider interval of energy up to 35 GeV the smeared experimental data have 
been recently obtained in Ref. [0. We represent these data and our NLO result in Fig. |3|, 
where we also show the parton model (PM) result as the dash-dotted curve. 

The D-function defined in the Euclidean region for positive momentum Q 2 is a smooth 
function and thus it is not necessary to apply any "smearing" procedure in order to have 
the possibility of comparing theoretical results with experimental data. We plot our results 
in Fig. f|, where we also show the experimental curve taken from Ref. || and parton model 
prediction. The shape of the infrared tail of the .D-function is sensitive to the value of the 
masses of the light quarks (the smeared Ra(s) function for A ~ 1-3 GeV 2 is less sensitive 
to the value of these masses). In our calculations we use the following masses m u = irtd = 
250 MeV, m s = 400 MeV, m c = 1.3 GeV, and = 4.7 GeV, which are close to the constituent 
quark masses and incorporate, therefore, some nonperturbative effects. Practically the same 
values of the quark masses were used to describe the experimental data in Refs. |T7|j25|] . 

4. In this paper we have applied the nonperturbative method of the a-expansion to 
describe the single-argument functions which are directly connected with the experimental 
data describing e + e~ annihilation into hadrons. An important feature of this approach is the 
fact that for sufficiently small value of the coupling (the ultraviolet region of momentum) it 
automatically reproduces the conventional perturbative results. We emphasize that the entire 
high-energy physics regime is accessible within this approach. Even going into the infrared 
region of small momentum, where the running coupling becomes large and the standard 
perturbative expansion fails, the a-expansion parameter remains small and we do not find 
ourselves outside the region of validity of the approach. 

We have considered two quantities that are convenient both for the theoretical analysis 
and for the model-independent comparison with experimental data. The first one is the 
Poggio-Quinn- Weinberg smeared function, Ra(s), defined in the Minkowskian region. The 
second one is the Adler .D-function defined in the Euclidean region. It should be emphasized 
that the smeared quantity Ra{<1 2 ) and the D-function have different sensitivity to QCD 
parameters. For example, in contrast to the smeared function taken with a reasonable value 



3 The PMS results correspond to a scale parameter A^T = 281 MeV that is rather small as 

compared with more modern value A^g ~ 370 MeV j24 |, which is compatible with R T - ratio used 
here. 
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of the parameter A the shape of the D-function in the infrared region is very sensitive to the 
values of the masses of the light quarks. At the same time, for A ~ 2-3 GeV 2 , the function 
Ra(s) is sensitive to the mass of charmed quark. Thus, these functions can test different 
parameters of the theory and, in some sense, are complementary to each other. 

It should be stressed that in this approach we use the parameters which are included 
in the Lagrangian and do not introduce any additional model parameters. Nevertheless we 
found good agreement between our results and the experimental data down to the lowest 
energy scale. At the same time, note that a straightforward attempt to apply the conven- 
tional operator product expansion, which leads to the "condensate" term in the D-function, 
a (a s G 2 )/Q 4 contribution, is not satisfactory because of the ill-definition of this expansion at 
small momentum. 
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FIG. 1. The smeared function Ra(s) for A = 1 GeV 2 . The solid line corresponds to the VPT 
next-to-leading order (NLO) result. The PMS result (the dashed curve) and the smeared experi- 
mental curve (the dotted line) taken from Ref. ^j. 
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FIG. 2. The smeared quantity Ra(s) for A = 3 GeV 2 . The solid, dashed and dotted curves are 
defined as in Fig. [I| The dash-dotted curve corresponds to the leading VPT order. 
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FIG. 3. The smeared quantity Ra(s) for A = 3 GeV 2 . The solid curve is the VPT result. The 
parton model (PM) prediction is presented by the dash-dotted line and the experimental curve from 
Ref. j^j is shown as the dotted line. 




FIG. 4. The plot of the D-function. The solid curve is the VPT result. The parton model (PM) 
prediction is presented by the dash-dotted line and the experimental curve from Ref. j^/ is shown as 
the dotted line. 
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